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ABSTRACT 

Various spectral emission lines from star-forming molecular cloud core L1517B manifest red asym- 
metric double-peaked profiles with stronger red peaks and weaker blue peaks, in contrast to the oft- 
observed blue-skewed molecular spectral line profiles with blue peaks stronger than red peaks. Invoking 
a spherically symmetric general polytropic hydrodynamic shock model for the envelope expansion with 
core collapse (EECC) phase, we show the radial flow velocity, mass density and temperature struc- 
tures of self-similar evolution for L1517B in a dynamically consistent manner. By prescribing simple 
radial profiles of abundance distribution for pertinent molecules, we perform molecular excitation and 
radiative transfer calculations using the publicly available RATRAN code set for the spherically sym- 
metric case. Emphatically, spectral profiles of line emissions from the same molecules but for different 
line transitions as well as spectra of closely pertinent isotopologues strongly constrain the self-similar 
hydrodynamics of a cloud core with prescribed abundances. Our computational results show that the 
EECC model reproduces molecular spectral line profiles in sensible agreement with observational data 
of Institut de Radioastronomic Millimetrique (IRAM), Five College Radio Astronomical Observatory 
(FCRAO) and Effelsberg 100 m telescopes for L1517B. We also report spatially resolved observations 
of optically thick line HCO + (l — 0) using the Purple Mountain Observatory (PMO) 13.7 m telescope at 
Delingha in China and the relevant fitting results. Hyperfine line structures of NH3 and N2H" 1 " transi- 
tions are also fitted to consistently reveal the dynamics of central core collapse. As a consistent model 
check, radial profiles of 1.2 mm and 850 /im dust continua observed by IRAM 30 m telescope and 
the Submillimeter Common- User Bolometer Array (SCUBA), respectively, are also fitted numerically 
using the same EECC model that produces the molecular line profiles. L1517B is likely undergoing 
an EECC shock phase. For future observational tests, we also predict several molecular line profiles 
with spatial distributions, radial profile of sub-millimeter continuum at wavelength 450/im, as well as 
the radial profiles of the column density and visual extinction for L1517B. 

Subject headings: dust, extinction — hydrodynamics — ISM: clouds — line: profiles — radio contin- 
uum: ISM — shock waves 



1. INTRODUCTION 

Observations over past two decades provide effec- 
tive diagnostics to examine dynamic properties of 
molecular clouds where star formations are taking 
place in relatively early phases. These observa- 
tions, including molecular spectral line profiles (e.g. 
Benson fc M^ersl [l989t ITafalla et all [200l lAguti et all 
2007f ) , radial profiles of (sub) mill i meter radio continua 
/e.g. iMotte Andre fe Neril [19981: IShirlev et all I2IM 
iStanke et all 120061: iNutter fc Ward- Thompson! [ 2007) , 
and i nfrared dust extinctions (e.g. lAlves. Lada fc Ladal 
HqoH), can be utilized separately and/or in combination 
to reveal several key hydrodynamic cloud features, such 
as core collapses, envelope expansions, travelling shocks 
and turbulence in the early evolution of star-forming 
cloud cores. 

On a relatively independent track, the theoret- 
ical model framework based on self-similar hydro- 
dynamic collapses has achie yed substantial develop- 
ment since late 1960s (e.g. iBodenheimer fe Sweigartl 
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1961 ILarsonl [T969I: IPenstonl [T969I: IShul [19771 IHunterl 



19771 ). and different variations of these theoretical 



mode l applications have been adva n ced in recent years 
/e.g. iTerebev Shu fe Cassenl Il98l ILou fe Shenl l200i 
iWang fe Loul 120081: iStahler fe Yen! I2009D . For grossly 
spherical molecular clouds, the gas flow velocity, mass 
density and thermal temperature structures of star- 
forming molecular cloud cores predicted by theoreti- 
cal models can be tested by several complementary 
high-resolution observations of selected source candi- 
dates. Attempts of this kind have been pursued all 
along by comparing theoretical cloud model results 
with the observed molecular spect r al lin e profiles (e.g. 
[ Zhou et all Il99l iTesti fc Sargentl Il99l ITafalla et all 
12004 120061 ) , (sub)millimeter r adio continuum emis- 
sions (e.s. IShirlev et all l200l IMotte fc Andrl 12001 
IHarvev. Wilner fc Mversri2003lL and infrared dust ex- 
tincti ons (e.g. lAlves. Lada fc Ladal 120011: iKandori et all 
1200511 

Molecular spectral line profiles are important to study 
physical structures and kinem atics for various phases i n 
star formation processes (e.g. IDvson fc Williams! 119971) . 
They need to be physically interpreted and can be prop- 
erly utilized to probe gas densities and temperatures 
through excitations of molecules, and detect cloud tur- 
bulence as well as systematic bulk flow motions through 
their line widths and Doppler shifts etc. Observations 
of spectral lines more frequently exhibit blue skewed 
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double-peak molecular profiles with stronger blue peaks 
and weaker red peaks (i.e. t he so- c alled blue profiles, see 
e.g. [Gregersen et all Il997l l2000t iLee. Myers fc Tafallal 
11999 ). These blue asymmetric spectral line profiles of 
molecular transitions are usually interpreted as signa- 
tures of core collapses with a static o uter envelope and 
a temperature variation involv e d (e.g. iZhou et alj|1993t 
IMvers. Evans fc Ohashil [200(1 IGao. Lou fe WuT l2009h' 
However, more and more high-resolution observations re- 
veal distinct red asymmetric double-peak profiles with 
weaker blue peaks and stronge r red peaks (i.e . the 
so-called red profiles, see e.g. iMardones et al.l IT9971: 
Thompson fc White! [2001 iFuller. Williams fe Sridharanl 



20051 IVelusamv et al.1 1200^ 7" which cannot be under- 
stood within the scenario of complete collapse models. 
The presence of 'red profiles' for molecular transitions 
most likely indicates a collapsing cent ral core surrounded 
by an expanding outer envelope (lLou fc Shenl 120041: 
Thompson fc W hite 20Ql IGao fc Loull2010ULou fc Gaol 



201 ID . though there have been alternative origins pro 



posed (see Gao & Lou 2010 and Lou & Gao 2011). 

The theoretical framework of self-similar hydrodynam- 
ics for envelope expa nsion with core co llapse (EECC) 
was first advanced by lLou fc Shenl ((2004) for an isother- 
mal gas cloud dynamics with spherical symmetry. This 
isothermal EECC model was later extended to a more 
general polytropic equation of state (EoS) with the spe- 
cific entro py conservation al ong streamlines (Wang & 
Lou 2007: rWang fc Lodl2008h . 

We advance the following scenario for the possible 
emergence of EECC dynamic phase in star-forming 
mole c ular clouds. Both theo retical (e.g. iBroderick et al.l 
2002t iStahler fc Yen! I2010D and observational (e.g. 



Lee fc Mversll2011[) works reveal that a static core would 



show expansive or oscillatory motions once being per- 
turbed on large scales. While a molecular cloud is under- 
going an expansive or oscillatory mode on large scales, its 
core may begin to collapse or contract due to nonlinear 
instabilities under self-gravity. With the enhancement of 
such central infalling dynamics induced by self-gravity, 
the collapsed or contracted region expands outward in a 
dynamic manner. Eventually, the entire cloud starts to 
collapse once the frontier of the collapsed region reaches 
the outer cloud edge. Thus, our EECC phase usually 
exists in the early stage of protostar formation, consis- 
tent with observations r egarding the intern al dynamic 
motions of starless cores (|Lee fc Myers 201lT) . 

Radiative transfer calculations based on such general 
polytropic EECC shoc k mod el have been systematically 
studied bv lGao fc Loul (f2010h and by Lou & Gao (2011) 
for the source FeSt 1-457 to emphatically demonstrate 
that the widely observed 'red profile' signatures in molec- 
ular line spectral profiles from low-mass star-forming 
clouds could be sensibly viewed as a diagnostic evi- 
dence revealing self-similar EECC hydrodynamic shock 
processes therein. To test the EECC interpretation of 
'red profiles', we investigate low- mass star- forming cloud 
cores with quasi-spherical configurations. In this paper, 
we focus on a candidate source of starless cloud core 
L1517B for comprehensive data fitting and comparison. 

To achieve high spatial resolutions, star-forming clouds 
with relatively simple internal structures yet without 
apparent envelope disturbances are chosen from the 
nearby Taurus complex. After a survey of known cases 



for star-forming cloud cores with red-skewed molecu- 
lar line profiles in the Taurus complex, we choose the 
cloud core L1517B for a more comprehensive investi- 
gation as it possibly involves global expansions. Cen- 
tral spectral line profile observations of molecular tran- 
sitions HCO+(l - 0), HCO+(3 - 2), H 2 CO(2i 2 - In), 
H 2 CO(2n - li ), CS(2 - 1) and CS(3 - 2) from this 
candidate source cle arly show red asymmetr ic spec- 
tral line profiles (e.g. iTafalla et all 120041 [20061) , which 
likely reveal the existence of an EECC shock dy- 
namic phase. Moreover, millimeter and sub-millimeter 
continuum mappings serve as independent constraints 
on r adial density and temperature profiles of L1517B 
(e.g. l Tafalla et al.ll2004HKirk. Ward-Thompson fc Andrei 
2005) for the same EECC shock model used for fitting 
molecular line profiles. 

Our motivation is to demonstrate specifically that the 
theoretical explanation of red profiles based on general 
polytropic EECC solutions w ith collapses, expansions 
and shocks (|Gao fc Loull2010l ). appears grossly consis- 
tent with available observations for starless cloud core 
L1517B. To further constrain and verify the EECC shock 
model, we present the result of a 5 x 5 (with 50" steps) 
spatially resolved HCO + (l —0) molecular line profile ob- 
servation using the 13.7 m telescope of Purple Mountain 
Observatory (PMO) at Delingha in Qinghai Province of 
China. We also report the model fitting results of these 
spatially resolved spectra on the basis of the same self- 
similar EECC shock model and show that the EECC 
shock dynamic phase will influence molecular line pro- 
files observed at different lines of sight (LOSs) from the 
cloud, rather than simply affects the central spectral line 
profiles alone. We also predict spatially resolved profiles 
for several other molecular emission lines, sub-millimeter 
continuum at 450/im, and the column density distribu- 
tion of L1517B for future observational tests. 

This paper is structured as follows. The background 
information and our motivation are introduced here in 
Section [TJ In Section [2l we specify model parameters 
and physical properties for L1517B using the framework 
of general polytropi c EECC hydrodynamic shock model 
([Wang fc L01JI20M) without magnetic fields. Model fit- 
tings of central spectra from several molecular line tran- 
sitions as well as spatially resolved spectral line profiles 
for the HCO + (l — 0) transition of L1517B are shown and 
analyzed in Section [3l We also venture to predict other 
spectral line profiles with spatial resolutions. In Sec- 
tion [H we show millimeter and sub-millimeter continua 
as further proof-tests to the same EECC shock dynamic 
model. Column number density and visual extinctions 
are predicted in Section 14.31 and we conclude in Section 
[5] with speculations. Details of analysis are summarized 
in Appendices A and B for a convenient reference. 

2. PROPERTIES OF STARLESS CLOUD CORE L1517B 

2.1. Earlier Studies on Cloud Core L1517B 

The candidate source L1517B as a star-forming cloud 
core is located in the Taurus c omplex, at an estimated 
distance of ~ 140 pc from us (|Eliad [l978h . According 
to the radio contin uum mappings and molec ular line in- 
tensity mappings (jTafalla et al.l I2002L l2006f ) , cloud core 
L1517B has a radius of ~ 150" (corresponding to ~ 
2.1 x 10 4 AU at a distance of ~ 140 pc) and a center iden- 
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tified at the right ascension aj2ooo = 04 h 55 m 18.8 s and 
the declination <5j2ooo = 30°38'04". Millimeter and sub- 
millimeter radio continuum data of L1517B have been 
acquired at 1.2 mm wavelength using the Institut de Ra- 
dioastronomie Mill i metri que (IRAM) 30m telescope by 
iTafalla et all ([2002I . I2004D . and at wavelengths of 850 /im 
and 450 /im using the Sub-millimeter Common-User 
Bolometer Array (SCU BA) on the James Clerk Maxwell 
Telescope (JCMT) by iKirk. Ward-Thompson fe Andrei 
( 2005). Model analyses of the 1.2 mm continuum have 
been done in two ways in the literature. One involves 
a presumed form of the number density radial profile 
no/[l + [ r / r o) a ] (where uq = 2.2 x 10 5 cm -3 is the cen- 
tral number density, r is the radius, tq = 35" is the 
reference radius and a = 2.5 is the scaling index) plus a 
constant du st temperature T = 9.5 K (e.g. Tafal la el all 
12004 120061) . The other is a static iso the rmal Bonnor- 
Ebert spher e (|Bonnorlll956tlEbertlll955T) by ITafalla et~aTI 
(|2002ll2004D . But the invoked Bonnor-Ebert spheres for 
clouds appear to be dynamically unstable given the in- 
ferred physical parameters from observations. 

These previous analyses offer relevant information for 
L1517B, with a constant temperature of ~ 9.5 K and a 
number density of ~ 2.2 x 10 5 cm~ 3 in the inner region. 
Information has also been drawn from the 850 /zm con- 
tinuum observation of L1517B (e.g. Kirk et al. 2005), in- 
cluding a temperature estimate of ^ 10 K, a full width at 
half maximum (FWHM) diameter of ~ 0.019 - 0.028 pc 
(by a 2D Gaussian fit to the 850 /im data), a central vol- 
ume number density of ~ 5 x 10 5 cm' 3 and a total cloud 
mass of ~ 1.7 M Q within the 150-arcsec aperture. Ra- 
dio continuum mappings on shorter wavelengths (850 /im 
and 450 /im in K irk. Ward- Thompson fc A ndre 20Qa) re- 
veal more substructures (thus not that spherically sym- 
metric as compared to the mapping of 1.2 mm continuum 
emissions) and may be used to constrain theoretical mod- 
els (see subsection 4.2). Polarization observations for 850 
/zm and 450 nm emission continua are performed usin g 
SCUBA by IKirk. Ward-Thompson fe Crutcherl (|2006l) . 
from which a magnetic field strength of ~ 30 /iG in this 
cloud core is inferred. For the moment, magnetic field is 
not included in our current model formulation. 

Various molecular transition lines from L1517B have 
been observed usi ng the IRAM 30 m telescope by 
ITafalla et all (|2004L 120061 ). and spectral line profile fit- 
tings based on a static gas sphere close to the center 
with an empirical density profile, a constant tempera- 
ture and an ad hoc outer (r > 10 4 AU) flow velocity gra- 
dient are performed therein. They also argued for other 
scenarios, such as rotation and pure contraction with cer- 
tain cooling mechanisms near the center to g enerate such 
asymmetric profiles (e.g. ITafalla et af]|2004D . There are 
inconsistencies in such empirical approach from the the- 
oretical perspective as also discussed in Gao, Lou & Wu 
(2009). Nevertheless, these model fittings might provid e 
gross information for L1517B. In ITafalla et all (|2004D . 
they declare the existence of internal motions of the or- 
der of ~ 0.1 km s _1 in addition to turbulence, which we 
shall see in the following analysis might be related to 
the cloud systematic infall and expansion motions pre- 
dicted by our general polytropic EECC shock hydrody- 
namic model. Another important result is the finding 
of molecular abundance drops towards the cloud center 



TABLE 1 

Six independent physical scaling parameters and 

DIMENSIONLESS SELF-SIMILAR EECC MODEL PARAMETERS 



Scaling h^' 1 


(km 


s -0.8) 


t (yrs) 


Value 


66.0 


3.0 x 10 5 


Parameters 


7 a 


A 


B b xi c 


Values 


1.2 


10.02 


4.58 1.87 



a 7 is the polytropic index of the EoS of a general polytropic gas. 

b Two integration constants A and B are mass and velocity pa- 
rameters in the asymptotic solution for x — y co [see eqns. JA11 1 in 
Appendix A] adopted as the 'boundary condition'. 

c xi is the dimcnsionlcss upstream location of an outgoing shock. 

(e.g. ITafalla et aLl l2006). which contributes to molecular 
line profile structures and will also be adopted in spectral 
line profile fittings of our model analysis (see subsection 
I3~4lfor details). 

2.2. Self-similar EECC Shock Model Parameters 

We adopt the general polyt ropic self-simi lar EECC 
hydrodynamic shock model of Wa ng fc Loul (|2008| ) to 
simultaneously fit both the (sub) millimeter radio con- 
tinuum emissions and several available molecular line 
profiles, especially those profiles with red skewed peaks 
and with spatial resolutions. The general polytropic 
model can semi-analytically and numerically describe 
self-similar behaviors of the hydrodynamic evolution, 
i.e., a time-evolving molecular cloud preserves its basic 
structural profiles (e.g. the density, temperature, veloc- 
ity and pressure profiles remain similar to their initial 
profiles), of a quasi-spherical polytropic gas under self- 
gravity with the specific entropy conserved along stream- 
lines (see Appendix A for basic nonlinear hydrodynamic 
partial differential equations (PDEs), the self-similar 
transformation, notations, and definitions therein) B A 
gas thermal temperature of ~ 10 K (e.g. I Tafalla et al.1 
120021 IKirk. Ward-Thompson fc Andrell2005l ) and a tvpi- 
cal cloud radius of ~ 0.019 - 0.028 pc Kirk et al. (2005) 
are considered for setting pertinent physical scalings in 
our self-similar EECC shock model. As L1517B appears 
in an early evolutionary phase of protostar formation, we 
may estimate a characteristic infall age of about ~ 3 x 10 5 
yrs for this cloud core according to Kirk et al. (2005) 
and iMverd (|2005f ) . We also adopt an outer radius of 
~ 2.1 x 10 4 AU for our cloud core model o f L1517B, by 
referr ing to the observed 150" radius (e.g. ITafalla et al. 
2002[) at an estimated distance of ~ 140 pc (e.g. lElias 

Using these empirical information of L1517B, we 
choose the self-similar model parameters and the phys- 
ical scalings of time t as well as sound parameter k as 
listed in Tables 1 and 2, of which the sound parameter k 
is derived from the following empirical length scale 

k l/2 t n g, 1Q 4 ATJ _ ^ 

The number density scale of L1517B cloud core for hy- 

5 To explore the origin of "red profiles" , we actually reduce the 
number of free model parameters in our data fittings by requiring a 
constant specific entropy everywhere (i.e. n + 7 = 2 is introduced), 
corresponding to the conventional polytropic cloud core. This con- 
straint is not necessary in the theoretical framework of our general 
polytropic model (Wang & Lou 2008). 
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TABLE 2 

Derived physical scaling k 2 and dimensionless self-similar 

EECC HYDRODYNAMIC SHOCK MODEL PARAMETERS 



Scaling 


k^ 2 (km s" 


-0.8) a 










(km 


Value 


66.8 














n 


q b mo c 


Vl 


ai 




V2 


a 2 d 




0.8 


0.039 - 


-0.17 


0.70 


1.85 


0.73 


1.55 





a The downstream sound parameter k 2 = (xi/x2) ki. 

b Paramctcr q = 2(n + 7 — 2)/(3n — 2) is determined by the poly- 
tropic index 7 and the scaling index n = 2 — 7 adopted in the 
self-similar EECC hydrodynamic shock model; this corresponds to 
a constant specific entropy everywhere. 

c The reduced central point mass mo is obtained by solving the 
self-similar hydrodynamic equations numerically and by matching 
the central free-fall solution (see Appendix A for details). 

d Parameters x± (regarded as independent in Table 1), X2, 01, 
V2, cti and a 2 are the dimensionless reduced upstream and down- 
stream shock locations, velocities and densities respectively, and 
are o b taine d from hydrodynamic shock jump conditions [sec cqns. 
||B2J| — ||B4|| ] in the self-similar shock model in Appendix B. 



drogen molecule H2 is taken as [eq 

' ' 10 4 cm 



U) in Appendix A] 
(2) 



where we adopt the mean molecular weight fj, = 2.29 as 
used by Harvey et al. (2003). This number density scale 
is one order of magnitude smaller than the empirical val- 
ues for the central num ber density ~2-4x 10 5 cm~ 3 
(e.g. lTafalla et al.ll2002t Kirk et al. 2005), as the density 
radial profile of our dynamic model increases rapidly to- 
wards the core center (see Fig. [T]). 

With these specified scalings, physical variables includ- 
ing radius, radial flow velocity, number density and ther- 
mal temperature of L1517B can be expressed in terms 
of the dimensionless self-similar vari able x and reduced 
dependent variables of x using eqns (| A5[i — (|A7[) . viz. 

(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 

The enclos ed m ass and the central mass accretion rate 
[see eqns. (|A6|) and (|A8|) in Appendix A] are 



n = 1.06 x 10 4 ir AU , 
r 2 = 1.07 x I0 4 x AU , 
U\ = 0.168 v(x) km s~ , 
ui = 0.170 v(x) km s 1 , 
Nx=N 2 = 3.49 x 10 3 a(a;) cm" 3 
Ti = 7.79 a(x)' y - 1 m(x) q K , 
T 2 = 7.98 a(x) 7 - 1 m(x) <z K . 



Mi 



°' 336 71 o\ M <3 

(3n — 2) 



and 



(10) 
(11) 
(12) 

where subscripts 1 and 2 refer to physical variables on 
the immediate upstream and downstream sides of the 
outgoing shock front, respectively. As required by the 
mass conservation across a shock surface, we confirm that 
Mi \ x=Xl = M 2 \x=x 2 = 2.35M by a direct numerical 



in \x ) 

M 2 = 0.345 — —r M Q , 
(3n — 2) 

M = 1.16 x 10~ 6 m Mgyr" 1 



i- 



Collapse region^^ 


Radial flow velocity 

R =1.98x10 4 AU ' 

R=1.68x10 3 AU 
In ' Shockwave 


Expansion region - — — ! 






Number density 


L1517B — -^^^ 

: M ,„r 3 - 89 M s > M =0.034 M s Shockwave 




Temperature 


Dynamic parameters 
" y=1.2, x,=1.87 
A=10.02, B=4.58 


shock wave 



R(10 AU) 

Fig. 1. — Structures of cloud core L1517B. From top to bottom 
are radial profiles of radial flow velocity in unit of km s _1 (positive 
values for radial infall), number density in unit of 10 4 cm -3 and 
temperature in unit of Kelvin (K), respectively. The abscissa is 
radius R in unit of 10 3 AU in a logarithmic scale. Infall radius 
R in {(= 1.68 X 10 3 AU) and shock radius -R sh (= 1.98 X 10 4 AU) ex- 
pand at speeds of ~ 0.021 km s — 1 and ~ 0.25 km s — 1 , respectively. 
Variables exhibit discontinuities across the shock front. Other pa- 
rameters for this general polytropic self-similar EECC shock solu- 
tion are summarized in Table 3. The radial density profile of our 
dynamic polytropic sphere with expanding envelope and free-fall 
collapsing core thus consists of two parts connected by an expand- 
ing stagnation surface, of which the inner part can be described as 
p <x r -3 / 2 [eq. IA12H 1 and the outer part as p oc r _2 / n [eq. (I A 1 11) 1 
(see the middle panel here). This broken pow er law radial densit y 
profile is a typical feature of starless cores (e.g. Casclli et al. 2002). 

computation with model parameters. To simplify physi- 
cal units of variables, we have already substituted valu es 
of fc, t and n as listed in Tables 1 and 2 into eq. (| A8|) . 

All dynamic model parameters in Table 1 are chosen 
in the procedure of data fitting of molecular spectral 
lines and dust emission continua from observations as 
described in Sections 3 and 4. In particular, red profiles, 
e.g. HCO+(l - 0) and HCO+(3 - 2) transitions, appear 
suggestive of an expanding envelope in L1517B; this mo- 
tivates us to invoke a self-similar polytropic EECC shock 
solution. Meanwhile, in order to fit the molecular spec- 
tral line profiles of these emissions, we introduce an out- 
going shock in the EECC model to avoid the outer en- 
velope of this molecular cloud expanding too fast. This 
shock emerges as an expanding flow rushes into the enve- 
lope (see the top panel of Fig. [T]) . We have systematically 
explored a wide range of model parameters (including 
polytropic index 7, asymptotic behavior characterized 
by mass and velocity parameters A and B, and location 
or speed of the outgoing shock) to identify the best-fit 
EECC shock model. The finally chosen EECC shock 
model with parameters in Table 1 does self-consistently 
and simultaneously fit (sub)millimeter emission continua 
and molecular spectral line profiles of L1517B. Velocity, 
density and temperature radial profiles at the present 
epoch of the star-forming L1517B described by the cho- 
sen polytropic EECC shock model are displayed in Fig. 
[TJ By the very dynamic nature, these radial profiles in 
Fig. [T]qualitat ively differ from those of b oth fitting mod- 
els adopted bv lTafalla etaLl (|2004 120061 ) and Kirk et al. 
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TABLE 3 

Physical properties of star-forming cloud core L1517B derived from our polytropic hydrodynamic EECC shock model 



Variables 


M 


iWtot a M {M Q yr" 1 ) b R inf (AU) R sh (AU) c «i (km s" 1 ) u 2 (km s' 1 ) d 


Values 


0.034 M© 


3.89 Mq 4.54 x 10~ 8 1.68 x 10 3 1.98 x 10 4 -0.03 0.12 



Note. — These physical pa rameters are grossly consistent with those estimated from SCUBA observations by 

(Kirk. Ward- Thompson & Andre 2005), except for the conspicuous feature of the variable temperature profile obtained self-consistently 
from our self-similar EECC hydrodynamic shock model, which differs from the static and isothermal assumptions in earlier models. 

(a) Mo and Mtot are the central point mass and the total cloud mass within a radius R = 2.1 X 10 4 AU. 

(b) Mo denotes the central mass accretion rate which decreases with time in our EECC shock model (see equation |A8t . 

(c) i?i n f is the expanding boundary separating the core collapse and envelope expansion regions and R s ^ stands for the shock radius. 

(d) u\ and u 2 are the upstream and downstream radial velocities across the outgoing shock front, with negative values for local inflows 



(2005). In principle, all these model fits should be con- 
strained simultaneously by other available observations. 

As fitted by the underlying hydrodynamic shock 
model, L1517B appears to involve a self-similar core col- 
lapse with an envelope expansion at a typical outflowing 
speed of w cxp ~ 0.1 km s _1 in the radial range of 1.7 x 10 3 
AU <i?<2.5xl0 4 AU, and a concurrent core collapse 
within the infall radius R ln f = 1.68 x 10 3 AU with a 
typical infall speed of U[ n f ~ 0.2 km s _1 . These infall 
and expansion speeds are in order-of-magnitude agree- 
m ent with the ~ 0.1 km s _1 internal speed estimated 
in iTafalla et al.l (|2004l ). There is an outgoing shock at 
Rsh = 1.98 x 10 4 AU with a travelling speed of ~ 0.25 km 
s . The central point mass which represents the mass of 
the pre-protostar is M ~ 0.034 Mq and the total cloud 
mass within a radius of 21000 AU is M tot ~ 3.89 Mq. 
Thus, there is not yet a star (in the sense of thermal 
nuclear burning) at the center of L1517B since the in- 
ferred mass of the pre-protostar is much smaller than 
the threshold mass value (~ 0.07 Mq) to initiate nuclear 
reaction. Besides, according to the standard equation for 
the accret ion luminosity L arr = GM+ M+ / R+ (e.g. equa- 
tion 2 in iKenyon fc Hartmannl 119951 ) , we can estimate 
an accretion luminosity of about ~ 0.02L Q for L1517B, 
which is likely below the curr ent infrared (IR) det ection 
limit. Moreover, according to iCaselli et~atl ([20021) . star- 
less cores are less massive (for example, a statistical mean 
mass is estimated as (M to t)— 3M Q , while the total mass 
of L1517B inferred from our model is ~ 3.89 Mq) than 
cloud cores with nuclear burning stars (e.g., (Mtot)— 
9Mq in contrast). Therefore, L1517B can be consis- 
tently considered starless, rather than a core with a star, 
as noted before. Our inferred total cloud mass is twice 
that derived from dust continuum flux densities by Kirk 
et al. (2005), bu t should still be classified as low- mass 
(M < 8 Mq, cf. IMcKee fc Osfaikerl [2001 star-fo rming 
cloud cores. The ratio of pre-protostar mass to the cloud 
mass is very low, at Mo/M to t — 0.87%, which implies 
that L1517B is in an early phase of protostar formation. 
The current central mass accretion rate of L1517B ob- 
tained from eq. (JT2J) is M = 4.54 x 10~ 8 Mq yr -1 , giv- 
ing an evolution timescale of t~E ~ Mq/Mq = 7.5 x 10 5 
yr. From the polytropic model with n < 1, the mass 
acc retion rate decreases with increasing time (see eq. 
IA8|) . which appears to be a t ypical feature in low-mass 
star formation sy stems (e.g. iSchmeia fc Klessenl 120041 : 
lEvans et al.1 120091 ) . Thus the realistic infall age of the 
cloud should be less than the derived evolution timescale 
t~E- Derived physical properties of the pre-protostellar 
core L1517B are summarized in Table 3 for a convenient 
reference. 



3. MOLECULAR SPECTRAL LINE PROFILES 

Molecular line profile observations reveal rich details 
for physical properties of protostellar cores, including 
important information of the cloud kinematics (e.g. 
iPavlvuchenkov et al.l [2008) . We perform spectral line 
radiative transfer (LRT) calculations for molecular line 
spectra us ing the publicly available num erical code 
RATRAN (iHogerheiide fc van der Takll2000l see website 
http://www.sron.rug.nl/ ~vdtak/ratran/ratran.html), 
which deals with both LK1' and non-local thermal 
equilibrium (non-LTE) excitations of molecular energy 
levels based on Monte Carlo method. This Monte Carlo 
code can handle both optically thin and thick lines, 
barring very large optical depths (e.g. > 100) caused by 
convergence problems. We apply the one-dimensional 
(ID) version of RATRAN code for spherical geometry. 
We adopt the EECC shock hydrodynamic model with 
parameters specified in Section [5] for molecular spectral 
line profile calculations. Here, the radial flow velocity 
u, number density N and gas temperature T are all 
derived self-consistently from this EECC model. The 
dust temperature Td in the cloud core is assumed to 
be equal to the gas temperature T as expected for 
high-density cloud cores where gas molecules and dusts 
have sufficiently frequent collisional excha nges (e.g. 
IGoldsmith fc Langerl 11971 IGoldsmithl l200l . When 
running the RATRAN code, we divide the spherical 
cloud core into 16 uneven shells with enough accuracy 
for calculations^ based on the principle that physical 
quantities in each shell do not vary significantly. We 
input the number density N, gas and dust temperature 
T, radial flow velocity u derived self-consistently from 
our EECC shock hydrodynamic model for radiative 
transfer calculations. Besides, the temperature of back- 
ground radiation is the cosmic microwave background 
(CMB) temperature T bg = 2.73 K. All molecular line 
transition data are obtained from th e Leiden Atomic an d 
Molecular Database (LAMDA, see lSchoier et al.l 120051) . 
and various molecular abundance ratios with respect to 
H2 molecules are presented in Table 4 and more details 
of abundance profiles can be found in subsection 13.41 

3.1. Red Skewed Double-Peak Line Profiles 

Numerical results (based on our EECC model) of cen- 
tral molecular spectral line profiles in comparison with 
relev ant observational data from ITafalla et al.l (|2004l 
2006) are displayed in Figs. [2] [3] and [4] Emission lines 

6 Actually, 12 and 14 uneven shells have also been tested sepa- 
rately in running the RATRAN code and one can sense the gradual 
convergence to the final results with little variance. 
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Fig. 2. — Central molecular spectral line profiles of L1517B 
with salient red skewed characteristics. The ordinate stands for 
the brightness temperature in Kelvin and the abscissa represents 
the LOS v elocity component. Histog rams are observational data 
taken from Tafalla ct al. (2004, 2006) while solid curves represent 
RATRAN LRT computational results based on our EECC shock 
model. Relevant abundances with respect to H2 are listed in Ta- 
ble 4. Our spectral profile fitting for the molecular transition line 
HCO + (l — 0) shows slightly stronger red peak than observed. How- 
ever, spatially resolved observation of this line transition reveals 
that the profiles in some spatial positions [see e.g. (—50", —50") 
in Fig. [5] exhibit somewhat stronger red peaks than our calculation 
in the central position due to non-spherical effects in L1517B. For 
molecular line transition CS(2 — 1), the intensity of the red peak is 
more or less the same bet ween the observatio n data and our fitting, 
while the fitting result of Tafalla ct al. (2004) is lower than the ob- 
servation data. These red asymmetric spectral line profiles reveal 
the presence of global expansions for the cloud core envelope. 



TABLE 4 

Rest-frame molecular line transition frequencies and 
molecular abundances in reference to h2 molecules 



Molecular Line 



Frequency 



Abundance 



Transitions 




(GHz) 






Xq 


HCO+(J= 1 - 


0) 


89.188523 


1.5 x 10 


-8 


(a) 


HCO+(J=3- 


2) 


267.557619 


1.5 x 10 


-9 


(a) 


H 2 CO(Jk-_ iK 


+1 = 2 12 - Lll) 


140.839502 


5.1 x 10 


-9 


(a) 


H 2 CO(J K _ lK 


+1 = 2n — lio) 


150.498334 


5.1 x 10 


-9 


(a) 


CS(J= 2-1) 


97.980953 


2.5 x 10 


-9 


(a) 


CS(J= 3-2) 




146.969026 


2.5 x 10 


-9 


(a) 


H 13 CO+(J= 1 


-0) 


86.754288 


7.5 x 10- 


11 


(a) 


DCO+(J= 3- 


2) 


216.112582 


4.2 x 10- 


10 


(a) 


SO(NJ= 23 - 


12) 


99.299890 


1.4 x 10 


-9 


(a) 


SO(NJ= 34 - 


23) 


138.178670 


1.4 x 10 


-9 


(a) 


N 2 H+(JFiF= 


101 - 021) 


93.176258 


8.0 x 10" 


10 


(b) 


N 2 H+(JFiF= 


344 - 233) 


279.511811 


8.0 x 10- 


10 


(b) 


NH 3 (JKFiF= 


1110.5 - 1111.5) 


23.694501 


2.5 x 10 


-8 


(c) 


NH 3 (JKFiF= 


2210.5 - 2211.5) 


23.722680 


2.5 x 10 


-8 


(c) 



Note. — Frequencies of molecular transition s are obtained from 
the LAMDA database (e.g. ISchoier et al.|[2005t) . 

a Molecules with central depletion hole. X = Xq for r > rj 10 i c , 
and X = 10~ 4 Xq for r < r^oic- F° r simplification and for keeping 
a minimum number of parameters in our model construction, we 
have assumed the same radius of abundance hole, r^oic = 5.1 X 10 3 
AU, for all molecular species in the fitting procedure. 

b Constant abundance without depletion hole. 

c The central abundance enhancement with Para-X(NHs) = 
Xo(NH3)[n(r)/no], where n(r) is the H 2 number density at ra- 
dius r. We adopt no = 2.2 X 10 5 cm -3 which is the same as 
Tafalla ct al. (2004) for comparisons. 



iTafalla et all ((20021 ) . we here take the observational data 
of ITafalla et al.l (I2004D in our model analysis^ Follow- 
ing ITafalla et all (|2004D . we produce the N 2 H+(1 - 0) 
and N2H + (3 — 2) spectra by averaging over a beam 
area of ~ 26" x 26" and ~ 11" x 11", respectively. 
As the NH 3 (J,iC) = (1,1) and (2,2) lines were ob- 
served simultaneously using the 100 m telescope of the 
Max Planck Institute for Radio Astronomy (MPIRA) at 
Effels berg near Bonn be tween 1998 October and 2001 
May (ITafalla et alJ 12001 ). we made an average over a 
~ 40" x 40" beam area to simulate relevant results. 
A constant intrinsic line broadening with Doppler b- 
parameter (defined as 1/e the half- width of a line profile 
in executing the RATRAN codeJH ~ 0.08 km s" 1 (i.e. 
with a FWHM of ~ 0.13 km s _1 for a Gaussian line pro- 
file) and a cloud receding velocity of u c i ou d ~ 5.73 km 
s _1 are used for all the data fitting procedure. 

Optical depths play consequential roles in producing 
asymmetric spectroscopic signatures in molecular line 



of HCO+(l - 0) and H 13 CO+(l - 0) are averaged over a 
beam area of about 50" x 50" to match with Five Col- 
lege Radio Astronomical Observatory (FC RAO) 13.7 m 
telesc ope observations in 2001 April (e.g. ITafalla et al.l 
120061 ). Spectral profiles for the following eight molecular 
transition lines, namely HCO+(3 - 2), H,CO(2i 2 - In), 
H 2 CO(2 u - 1 10 ), CS(2 - 1), CS(3 - 2), DCO+(3 - 2), 
SO (23 — 12) and SO (34 — 23), are averaged over a typical 
beam area of about 20" x 20" to match with the data 
acquired using the IRAM 30 m telescope in Spain be- 
tween 1999 Octo ber and 2002 November (|Tafalla et alJ 
120021 12004 I2006D . Noting that the observational data 
of the molecula r trans ition line CS(2 — 1) published 
in ITafalla et alJ ((2004) actually differ from those of 



7 In our fittings of spectral profiles from emission lines of CS, 
w e note that th e obse rved intensity of the red peak of CS(2 — 1) 
in ITafalla et all l |2002! ) appears too low to simultaneously fit the 
spectral profiles of CS(2 — 1) and CS(3 — 2) tr a nsitio n lines. We 
then checked the observations of Tafalla ct al. ( 2004) and found 
that Tafalla et al. (2004) have made a correction of the relevant 
data, and the new data appears more consistent with our model 
fitting calculations. The difference between the CS data of Tafalla 
et al. (2002) and Tafalla et al. (2004) may be caused by different 
radio telescopes between FCRAO 14m observations (Tafalla et al. 
2002) and IRAM 30m observations (Tafalla ct al. 2004). 

8 Both the thermal (At>x) and turbulent (Atit) components that 
contribute to the line broadening have been included in RATRAN 
code calculations. An-p is estimated by ~ {k^T /m) 1 / 2 with T 
consistently obtained from each shell and m is the mean parti- 
cle weight, and Aft = ex (Doppler b-parameter). The total line 
broadening Av = [(Av T )' 2 + (Av t ) 2 ] 1/2 . 



EECC for Star-Forming Cloud Core L1517B 



7 



profiles. Optically thick molecular line transitions of 
HCO + and CS display deeper self-absorption dips than 
optically thin emissions from molecular line transitions 
of H2CO, which have almost no central dips and show 
only stronger red shoulders (see Fig. [5]). Moreover, the 
HCO + (3 — 2) line transition manifests less apparent self- 
absorption dip in comparison with HCO + (l — 0) tran- 
sition. This contrast is likely due to lower energy level 
populations on J=3 and 3—2 in such a cold interstellar 
medium (ISM) environment of T ~ 10 K, leading to a 
lower optical depth and s ource function of the J= 3 — 2 
line transition in HCO+ (jGao fc Loull2010f ). 

Emissions from the same molecule but for different 
energy level transitions strongly constrain the under- 
lying polytropic hydrodynamic cloud core model, be- 
cause all the dynamic and thermal parameters as well 
as molecular abundance profiles should remain identi- 
cal in the LRT calculations using the RATRAN code. 
Fig. H shows the spectral profiles of HCO+, H 2 CO 
and CS each for two distinct transitions from central 
L1517B. These molecular line spectral profiles present 
explicit red skewed double-peak signatures (i.e. red 
profiles), consistent with the plausible existence of self- 
similar E ECC shock phase (Lou fc Shen 2004; Shen & 
Lou 2004: lThomnson fc Whitell200i iGaofc LouII201QD F1 
Such a cloud envelope expansion with core collapse may 
arise from either l arge-scale cloud radial oscillations (e.g. 
lAguti et al.ll2007[ ) or molecular out flows driven by pro- 
tosta r embedded in the core (e.g. Thompson & White 
2004). It has been claimed that the presence of molec- 
ular outflows is always associated w ith the evidence 
of non-Gaussian CO line wings (e.g. iThompson et al.l 
12001 . However, spectra for C 18 0(1 - 0), C 18 0(2 - 1), 
C 17 O(l-0) and C 17 Q(2 - 1) towards L1517B (see figure 
8 of iTafalla et al.l [20041 ) do not manifest such evidence 
of outflows. Besides, no spatially separated "molecular" 
outflow lobes has yet been detected in this star forming 
cloud core. We thus propose that the global envelope 
expansion with central core collapse revealed in L1517B 
may originate from damped acoustic radial pulsations 
on large spatial and temp oral scales in mol e cular cloud 
cores (Lou fc Shen 2004; lLada etTal] 120031: iKeto et al.l 
f2006t IGao fc LoullMfjl: ILou fc Gadl20H Fl 

3.2. Other Relevant Molecular Transition Lines 

As consistent checks and further constraints, we 
also compute central spectra of isotopologues, namely 
H 13 CO+ and DCO+, of the formyl cation HCO+ and 
two distinct line transitions of sulfur monoxide SO using 
the same EEC C shock model in co mparison with obser- 
vation data of Tafalla e t al.l (|2006l ). Fittings of isotopo- 
logue emissions with constant abundance ratios to each 
othei[3 provide an effective way to investigate the in- 
fluence of abundance patterns on spectral profiles. As 
already noted in subsection 13. H emissions from identi- 
cal molecules with different energy level transitions serve 
as effective evidences to verify our EECC shock model. 

9 A schematic explanation for the origin of red pr ofiles in molec- 
ular cloud cores can be found in Lou & Gac] pOlU ). 

10 More explanations for the coexistence of expansion and col- 
lapse are elaborated in introduction and discussion sections. 

11 We simply multiply a constant ratio to the molecular abun- 
dance of HCO^" for its isotopologues in our model profile fittings. 



In extensive numerical explorations, we adopt molecular 
abundances and isotopic r atios which a r e som ewhat dif- 
ferent from those used in ITafalla et al.l ()2006l ) . e.g. the 
SO abundance with respe ct to hydrogen mole cule H2 dif- 
fers from that adopted by Taf alla et al.l (|2006| ) due to the 
variation in the chosen radius of central depletion hole. 

There is no conspicuous appearance of red skewed 
double-peak profiles in molecular line emissions in Fig. 
[3j instead, single-peak lines are observed in these opti- 
cally thin lines as expected. As discussed in subsection 
13. H the absence of red asymmetries and self- absorption 
dips in spectral line profiles is attributed to the decrease 
of molecular level populations (which in turn leads to 
the optical depth decrease). This point is highlighted 
by comparing the spectral line profile of the optically 
thin H 13 CO+(l - 0) in Fig. [3] with that of the opti- 
cally thick HCO+(l - 0) in Fig. [2 [n.b. the optical 
depth of H 13 CO+(1-0) is ~ 1/20 that of HCO+(1-0)]. 
Therefore, the combined examination of optically thick 
and thin molecular lines offers an important diagnosis to 
constrain la rge-scale structures of low-mass star forming 
cloud cores (IGao fe Lodl20lol: ILou fc Gadl2lHTlh 

Compared with observations, our model results in Fig. 
02 especially those of DCO + and SO, tend to be broader 
in line profile widths. Actually, model lines of these 
molecules could be narrowed by reducing their depletion 
hole size43. For illustrations, we have decreased the de- 
pletion hole radius rhoic from 5.1 x 10 3 AU to 2.2 x 10 3 AU 
for both DCO + and SO. Meanwhile, their abundances 
are reduced to 1.5xl0~ 10 and5xl0~ 10 , respectively, and 
then both molecular lines become narrower (see dashed 
curves in Fig. |3|). Physical explanations for the reduc- 
tions of depletion holes are as follows: The difference be- 
tween the abundance hole radius of the DCO + and that 
of the HCO + results from a central increase in the deu- 
terium fractionation caused by the CO depletion, which 
partly co mpensates the DCO + freeze out at the inner 
core fe.g. ITafalla et^l200fih . The non- C arbon-bearing 
molecule SO does not suffer from the depletion of Carbon 
components, so that its depletion hole may differ from 
those of Carbon-bearing molecules. Here, SO provides 
complementary information to emissions from Carbon- 
bearing molecules for the protostar-forming cloud core. 

3.3. Molecules Without Central Depletion Holes 

Molecules with central depletion holes offer important 
diagnostics to probe the dynamic structure of the outer 
molecular envelope of L1517B. Meanwhile, molecular 
tracers of the central core, such as N2H 4 " and NH3, offer 
essential information to examine the collapsing core, as 
these molecules do not freeze out onto dust grains at typ- 
ical core densities and th ey are present in th e gas phase 
throughout the core fe.g. ITafalla et aLll2004|) . As shown 
in Table 4, a constant N2H 4 " abundance and a variable 
abundance profile with a central enhancement for NH3 
have been adopted in our RATRAN calculations. 

12 The reduction of the depletion hole (still larger than the col- 
lapsing core) would bring extra emission contributions from the 
denser and hotter inner region, so we need to decrease the total 
abundance of these molecular species to better fit observations. 
Due to the lower abundance, emissions from the outer expanding 
envelope which mainly contribute to the edges of the line spectra 
are weakened, and this leads to narrower spectral line profiles. 
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Fig. 3. — Molecular line profiles of central spectra for four transi- 
tions H 13 CO+(1-0), DCO+(3-2), SO(23- 12) and SO(34-23) . 
Histograms are observational data taken from[Tafalla ct al] l|2006fl 
while solid curves represent LRT fitting calculations based on 
the same self-similar EECC shock model. Solid curves represent 
the model results with the universal depletion hole size rhole = 
5.1 x 10 3 AU for all species. Constant isotopic ratios of 12 C/ C 
and H/D (see Table 4 for details) are adopted in our model fitting 
analysis. Dashed curves displayed in the DCO+ and SO line pro- 
files correspond to the simulation results for a smaller depletion 
hole with r^oic = 2.2 X 10 3 AU, and lower molecular abundances 
as described in the text. The absence of red asymmetry and self- 
absorption dip in these molecular spectral line profiles is attributed 
to the decrease of optical depths because of the lower energy level 
populations for these molecular line transitions. 
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Fig. 4. — Molecular line profiles from the central part of L1517B 
without depletion holes. Multi-peaks present in both molecular 
transitions are caused by hyperfine (hf) splittings rather than by 
the cloud co re dynamic co ll apse. Histog rams are observational data 
taken from ITafalla et all ||2002, 2004). As shown in Table 4, a 
constant N2H+ (top panels) abundance and a central abundance 
enhancement pattern for NH3 (bottom panels) have been adopted 
in order to fit the observational data. 



Unlike other molecular transitions discussed in this pa- 
per, rotational levels of N2H 4 " are split into multiple 
hype rfine (hf) components by the two nitr ogen atoms 
(e.g. ICaselli et alJ 12001 ITafalla et all 12001 . Such hf 
splittings make the radiative transfer computation much 
more complicated. Fortunately, available N2H" 1 " molec- 
ular data file containing hf structures is provided in 
LAMDA, which simplifies our calculations considerably. 
In all, 15 different transitions of N2H + (with 7 distinct 
frequencies due to the overlap of some transitions) for 
,7=1 — and 29 different transitions (with 26 distinct 
frequencies) for J = 3 — 2 are actually involved in their 
spectra. We first compute populations of these hf sub- 
levels and related radiative transfer processes with regard 
to each hf transition, and then superpose all the individ- 
ual hf lines to compare with observations (see Fig. 2]). 
No obvious line broadening deviation from the observed 
lines has been found in our numerical fittings (either in 
J = 1 — or J = 3 — 2). While a significant contract- 
ing central region with inward motions of the order of 
~ 0.2 km s _1 has been predicted by our EECC model, 
this region is very small as compared to the entire col- 
lapsing core (~ 1/27 in volume) and is almost negligible 
compared to the entire cloud (~ 1/27000 in volume). 
Therefore, this region contributes little to the broaden- 
ing of molecular line profiles. Such a fairly "high" speed 
of collapsing region is possible in a cloud since gas mate- 
rials immediately above the surface of the central dense 
core (or pre-protostar) approach a state of free-fall. 

The ammonia molecule consists of the so-called or- 
tho and para species, which coexist almost indepen- 
dently of each other because normal radiative and col- 
lisional transitions do not change spin orientations (e.g. 
IHo fc Townei[l98l . The observed NH 3 (J,iT) = (1,1) 
and (2, 2) inversi on lines arise from para-NrTa (e.g. 
ITafalla et al.l [20 02) . so we only focus on this particular 
species in our modelling. Like N2H+, ammonia NH3 also 
has hf splittings due to the interaction between the elec- 
trical quadrupole moment of the nitrogen nucleus and the 
electric field of electrons. Apart from this major con- 
tribution, other weaker interactions, including the In - 
J (where In and J are the nitrogen spin and the to- 
tal angular momentum of ammonia, respectively) and 
I- J (where I is the sum of hydrogen spins) magnetic 
interactions, as well as H-N and H-H spin-spin inter- 
actions, further split th e transition compon ents on the 
order of - 40kHz (e.g. IHo fc Townei [1981 . Unfortu- 
nately, molecular data with hf splittings for radiative 
transfer calculations a re not availab l e in L AMDA. We 
follow the approach of Taf alla et al.l ([2002D by separat- 
ing the calculation into two steps, namely, the solution 
of the level excitation and the prediction of the emergent 
spectrum. Using the RATRAN code, we first adopt the 
ammonia data file without hf splittings from LAMDA to 
obtain the combined population of each rotational en- 
ergy level. We then compute the population of each 
hf subleveJ3 by assuming that the sublevels are popu- 
lated according to their statistical weights. Once pop- 
ulations for the hf sublevels are derived, we take the 

13 Information about hf sub-levels are taken from the NASA-JPL 
molecular database. Only the electrical quadrupole hf energy level 
structure has been considered in the data files. Therefore, 5 distinct 
components of the transition frequency for both (J, K) = (1, 1) and 
(2, 2) are taken into account in our radiative transfer modelling. 
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full hf structures into account and predict the emergent 
(J,K) = (1,1) and (2,2) spectra (as shown in Fig. 0) 
by integrating the radiative transfer equation along the 
LOS and summing up each hf componcnt P"^. The validit y 
of this approach was discussed by ITafalla etaLl (|2002D . 
We briefly comment here: according to quantum statis- 
tics, the population ratio between two different levels is 
given by N 2 /Ni = (32/51) x e^ AE ^ kBT \ where gi are 
the statistical weights (degeneracies) for level i (i=l, 2) 
and AE — E2 — E\ is the energy difference between the 
two energy levels. As long as AE <C k B T, i.e. the energy 
difference is far less than the energy scale involving the 
excitation temperature (which is usually satisfied when 
we consider the population of hf sub-levels for molecules 
in starless cores), we can then approximately take the 
population ratio as simply 32/31- The excitation tem- 
perature and the kinetic temperature can be equal under 
thermodynamic equilibrium, but this is not necessarily so 
in general. 

Ammonia NH3 spectral data are usually used to 
probe the central gas kinetic temperature of star-forming 
molecular cores (e.g. ITafalla etaLll200l . The database 
of dense molecular cores mapped in the (J,K) — (1,1) 
and (2, 2) transition lines of NH 3 was co mpiled and pre- 
sented by IJiiina. Myers fe Adamsl (fl999), who conclude 
that the temperature distribution of molecular cores in 
Taurus complex is very narrow around ~ 10 K. In our 
polytropic EECC shock model, however, the gas temper- 
ature gradually rises up towards the center and reaches 
~ 15 K in the collapsing core of L1517B as shown in 
Fig. [Q According to Tafalla et al. (2004) and Ho & 
Townes (1983), the gas temperature of ~ 10 K inferred 
from NH3 observations relies on the presumed (static) 
density profile, which is an empirical asymptotic power- 
law envelope with a "flat" central region. It appears 
that such a constant ~ 10 K temperature reasonably 
fits the gas kinetic temperature for L1517B derived from 
NH 3 data analysis (e.g. figure 4 of ITafalla etaLl 120041 ) . 
We note, however, that the gas kinetic temperature thus 
determined is inferred dir ectly from the rotational tem- 
perature Tl 1 of NH 3 (e.g. IWalmslev fe UngerechtsHl983t 
ITafalla et al.l 120041 ) which relates to the relative bright- 
ness temperature/population betwe en the (J, K) = (2, 2 ) 
and (1, 1) radiative transitions (e.g. lHo fc TownesH l983): 
therefore, the gas kinetic temperature thus derived only 
represents an "averaged" temperature along the LOS, 
not necessarily indicating the actual temperature at a 
certain spatial point. In other words, it is possible that a 
radially variable temperature profile could also reproduce 
the NH3 data as shown by our dynamic model fitting 
analysis (see also Galli et al. 2002 for variable tempera- 
ture profiles in molecular cloud cores). Noting that the 
temperature profile of our model reaches ~ 15 K o nly 
in the very inner region (~ 1000 AU from the centciFT) 
while remaining around ~ 10 K over a wide radial range 
in the outer portion (see Fig. [1]), we would expect that 
our radially variable temperature profile may mimic the 

14 We compute the Einstein A c oefficients for differe nt hf line 
transitions based on equation (9) of Pickett et al. (1993). 

15 The gas kinetic temperature for L1517B deriv ed from NH3 
obser vations starts from ~ 900 AU (see figure 4 of ITafalla et al.l 
2004); therefore, the "hot" region in our model occupies the inner 
core and contributes only a very small part to the LOS averaged 
gas kinetic temperature. 



constant LOS temperature distribution of the NH3 anal- 
ysis for L1517B. Besides, as the temperature d istribution 
of L1517B in figure 4 of ITafalla etaLl (|2004D represents 
LOS averages spatially smoothed with 40", these aver- 
ages over the telescope beam resolution would further 
flatten the temperature distribution and make the whole 
profile closer to a roughly constant ~ 10 K distribution. 
Moreover, as our EECC model shown in Fig. U as a 
whole reasonably reproduces both numerous molecular 
line emissions (including those of NH3 and N2H" 1 ") and 
(sub)millimctcr dust continuum observations of L1517B 
(see Table 4 and Section 4 for details), we conclude that 
this hydrodynamic self-similar EECC model is a viable 
physical description of cloud core L1517B. Our temper- 
ature variation would only introduce minor corrections, 
because the inner "hot" region is fairly small compared 
with the entire molecular cloud (i?i n f < 0.1i? ou t)- In 
other words, a variable temperature profile might ap- 
pear somewhat flattened in this type of inferences by 
LOS NH3 observations. Noting that our radial temper- 
ature profile drops below ~ 10 K in the outer envelope 
portion, it is possible to produce a 10 K "averaged" rota- 
tional temperature from our EECC model. Jijina et al. 
(1999) have adopted the T^-Tk (here Tk is the gas ki 



netic temperature) from I Walmslev fc Ungerechtsl ( 1983f ) 
to convert the inferred rotational temperature into gas 
temperature, which may cause further uncertainty in the 
determi nation of cloud temperature sin ce the T^-Tk re- 
lation of lWalmslev fe Ungerechtsl ()1983f ) may not be that 
accurate (e.g. ITafalla et aLll2004l ). 

3.4. Influence of Molecular Abundance Distributions 

The spatial distribution of molecular abundances plays 
a crucial role in producing molecular spectral line profiles 
and thus h as received considerable at t ention in the liter- 
ature (e.g. iHerbst fc Klempered [19731 : iRawlings fc Yatesl 
I200H iTsamis et al.ll2008D . Also" the depletion of molecu- 
lar species by adhesion onto cold dust grain surfaces has 
been discov ered in centra l regio n s of s t ar- forming cloud 
cores (e.g. ITafalla et all 120021 l2004t IWalmslev et~aTl 
l2004f ). We have selected molecular abundance ratios (rel- 
ative to the number o f H2 molecules) by referring to 
ITafalla et all (|2004 120061) yet with variations to probe 
molecular distributions in L1517B. Central molecular de- 
pletion holes with very much lower abundance ratios 
[i. e. ~ 10~ 4 o f the ratio in the outer cloud layers as 
in ITafalla et al.l (e.g. I2006f )] are essential in reproducing 
absolute intensities and relative strengths of 'blue' and 
'red' peaks in molecular spectral line profiles. It is pos- 
sible to treat the radii of abundance holes, which can 
vary for different molecular species in general, a s extra 
independent parameters as in lTafalla et al.l (120061 ). How- 
ever, in order to focus on exploring the dynamic process 
of L1517B and to reduce the number of independent pa- 
rameters, we have assumed a common size of depletion 
hole for all molecules, except for NH3 and N2H + molecu- 
lar transitions. Our chosen abundance distributions and 
the radius of the depletion hole rhoic used in the data 
fitting are summarized in Table 4 for reference. Based 
on LRT RATRAN calculations, intensities of both 'blue' 
and 'red' peaks would decrease with the increase of Thoic! 
this is also intuitively sensible. In fact, the intensity of 
the stronger 'red' peak appears to be much more sensi- 
tive to the variation of rhoie than that of the weaker 'blue' 
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Fig. 5. — We show here the grid-map for the observed molecular transition HCO + (l — 0) spectra with a 50" step for L1517B. The emission 
intensities are expressed in main beam temperature T m b (using Kelvin as the unit) with elevation corrections f(6) = 1 — 0.06(±0.02) cos 8 
where 9 is the elevation angle of the Delingha telescope, and the efficiency of the Delingha 13.7 m telescope is ?7 m b = 0.62. This 5x5 grip 
map is centere d at 0,12000 = 04 h 55 m 18.8 s , 5j2000 = 30°38'04" based on the central position of the 1.2 mm continuum source of L1517B 
as observed by[Tafalla ct al. (200!|). The LOS velocity resolution is Av = 0.04 km s — 1 , and the root-mean-square error in the antenna 
temperature T£ is <r rms = 0.083 K. The legend panel to the upper right corner illustrates the central panel as an example. The spectral 
profiles are not perfectly spherically symmetric with the eastern spectra slightly shifting towards the blue. The difference between the 
eastern and western spectra might be due to the rotation of L1517B about its north-south axis. However, double-peaked red profiles are 
evident in HCO+(l - 0) lines across the cloud globule and an average over the same distance from the center reveals stronger red peaks in 
all radial positions (see Fig. |S}. In this sense, the presence of red profiles is indicative of cloud core expansive motions (see also lAguti et al.1 
120071) . 



peak does. This may allow us to estimate the radius of 
the depletion hole in L1517B by carefully calibrating the 
'blue' to 'red' peak intensity ratios in pertinent molecular 
spectral line profiles. 

3.5. Spatially Resolved Spectral Line Profiles 

Central spectral profile model fittings of various molec- 
ular line transitions provide an important and effective 
approach to probe both dynamical and thermal struc- 
tures of pre-protostellar cores. However, only presenting 
signatures from the cloud center, central spectral pro- 
file fittings of molecular lines might still have certain 
ambiguity and uncertainty in determining all pertinent 
physical properties of star-forming cloud cores. Mapping 
molecular transition line profiles with high enough spa- 
tial resolutions, on the other hand, would further con- 
strain the large-sc ale physical condi t ions within molec- 
ular clouds (e g |AgutielaL| 120071: iGao fe Loul 120101: 
lLou fc Gaol 120 111 ). Here, we report a position-switch 
observation (5x5 grids) of spectral profiles for the 



HCO + (l — 0) line transition, which is the most prominent 
diagnostics for dynamic motions within L1517B, using 
the 13.7 m telescope at Delingha in Qinghai Province of 
China and present our model fitting results of spatially 
resolved spectral profiles based on the same underlying 
EECC shock model. Meanwhile, we predict other perti- 
nent red skewed spectral profiles of molecular lines with 
spatial resolutions for future observational tests. 

3.5.1. Data Acquisition and Comparison 

We recently observed L1517B for HCO+(l - 0) at 
89.186769 GHz using the 13.7 m millimeter- wave radio 
telescope of PMO at Delingha during April 12-22, 2010. 
The telescope is ~ 3200 m above sea level with an ex- 
tremely dry air. One SIS receiver operating at ~ 90 GHz 
was used in the observation together with a FFTS digital 
spectrometer with 16384 channels and a working band- 
width of 200 MHz, which gives a velocity resolution of 
0.04 km s^ 1 at the frequency of 89 GHz. The position- 
switch mode was chosen, with the pointing and track- 
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Fig. 6. — A direct comparison between the central spectra, which 
have been generated from the average of typically 5 spectra within 
a 50" radius from the core center, of the molecular transition 
HCO + (l— 0) from L1517B observed by Delingha 13.7 m millimeter- 
wave radio telescope in mid-April 2010 (solid histogram) and by 
FCRAO 13.7 m telescope in April 2001 (dashed histogram). The 
ordinate stands for the main beam brightness temperature T m t, 
converted using the beam efficiencies (0.62 for Delingha telescope 
and 0.55 for FCRAO telescope respectively) from the respective an- 
tenna temperatures Tf . The Delingha observation has a velocity 
resolution of 0.04 km s^ 1 as compared with the velocity resolution 
between 0.03 km s -1 and 0.07 km s — 1 of the FCRAO observation, 
which used the QUARRY array receiver in frequency switching 
mode together with the facility correlator. With longer integration 
times for the observation of each spatially resolved regions, the re- 
sult from the Delingha observation has relatively higher signal to 
noise ratio (S/N) compared with the FCRAO data. 

ing accuracy within the range of 2" — 7". The telescope 
beam size is approximately 50" at this frequency, and we 
adopted a mapping step size of 50" in observation. The 
typical system temperature T sys is ~ 250 K during our 
observations. The standard chopper wheel calibration 
was used during the observation runs to get the antenna 
temperature T£, which has been corrected for the atmo- 
spheric absorption and telescope elevations. The beam 
efficiency of 0.62 is used to convert the telescope antenna 
temperatures Tj£ into main beam brightness tempera- 
tures T m b. The net integration time for most positions 
is 60 minutes, and is 80 minutes for the central position, 
resulting a root-mean-square (rms) noise of cr rms = 0.083 
K for the measured antenna temperature T£. 

The spectral profile data are processed with the stan- 
dard CLASS software. Linear baselines are removed 
from all spectra and an elevation correction f{9) = 
1 — 0.06(±0.02) cos (9 (9 is the elevation angle) is adjusted. 

We also made an average of five spectra within a 
50" radius from the cloud core center (i.e. the central 
spectrum together with four spectra that are 50" away 
from the cen t er) to directly compare with the result of 
iTafalla et all ()2006[ ) obtained by the FCRAO 13.7 m tele- 
scope in April 2001. From the data comparison of Fig. 
El we see that these two observational results are gen- 
erally consistent with each other with the notable ex- 
ception of integrated intensities around the dip between 
the two peaks. This difference may be tolerated within 
uncertainties of the observations and might be due to 
slight discrepancies in telescope calibrations or velocity 
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Fig. 7.— The HCO+(l - 0) grid -map of four spectral profiles 
at a distance of ~ 20" around the central spectral profile plot for 
L1517B (see Fig. [S]for the average of these four spectra). Emission 
intensities are expressed in the main beam brightness temperature 
T m b- The net integration time for each spectral profile with the 
impact parameter b= 20" is 50 minutes, and is 80 minutes for 
the spectral profile at the center. While these 20" spacing spectra 
overlap previous observations with the beam size of 50", it is still 
valuable to examine global dynamic properties of L1517B from 
these spectra, as the effect of possible rotation about the north- 
south axis in the inner cloud region appears weaker than in the 
outer region. Compared with Fig. \5\ apparent red profiles in nearly 
all these inner spectra (n.b. the red asymmetry is not particularly 
obvious in the eastern spectrum) strongly indicate that this red 
asymmetry is an intrinsic characteristic of L1517B caused by its 
internal motion rather than by rotation. 

resolutions. This comparison of data from two indepen- 
dent observations 9 years apart confirm the validity of 
both observational data from FCRAO 2001 April and 
Delingha 2010 April. 

To examine whether these "red-profiles" are caused by 
rotation, we have observed four spectral profiles at a dis- 
tance of 20" from the center (Fig. [7]) since the rotational 
effect is much weaker in the inner region than in the outer 
region. Conspicuous red asymmetric profile characteris- 
tics in nearly all these spectra exclude the possibility that 
a pure rotation produces these "red-profiles" . 

3.5.2. Model Fittings of Observational Data 

Simultaneous comparisons between model results and 
observational data for molecular line profiles with dis- 
tinct impact parameters b (distance of LOS from the 
core center) is another powerful way to examine the 
overall dynamical and thermal structures of star-forming 
cloud cores. As our model is spherically symmetric while 
the spatially resolved observation (Fig. [SJ exhibits non- 
spherical effects, we make an average over grids of data 
with the same radial distance from the center. We carried 
out observations of four spectra spaced by 20" from the 
core center under the same observational conditions (yet 
with 50 minute integration time each) to check EECC 
shock model (Fig. [7]). All the dynamic parameters as well 
as the abundance distribution and the intrinsic broaden- 
ing in these simulations are consistent with those used in 
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fitting the central spectrum of HCO + (l — 0) observed by 
the FCRAO 13.7m telescope (see Section 2). 

From the model fitting results shown in Fig. [8j we see 
that our computed molecular line profiles based on the 
general polytropic EECC shock hydrodynamic model co- 
incide with the observational data. In comparison with 
observations, molecular line profiles produced from our 
numerical model calculations show somewhat stronger 
red peaks and weaker blue peaks at the very center (Fig. 
[8j. This likely originates from the spatial asymmetry of 
the spectral line profiles, which might be caused by the 
possible rotation of L1517B about its north-south axis. 
The small blue shifts in the eastern spectra (see Fig. [5]) 
will reduce the intensity of the red peaks and enhance the 
blue peaks when we make an average of the observational 
data over the same radial distance. Meanwhile, these av- 
eraged spectra exhibit increasing broadening towards the 
outer envelope while the broadening of our numerical re- 
sults remain invariant. This phenomenon might be possi- 
bly attributed to the incremental influence of turbulence 
with increa sing radius, which is typical in star-forming 
clouds (e.g. IMcKee fc T an 2003). In order to reduce the 
number of free parameters and to better understand the 
dynamics of L1517B, we simply assumed that the intrin- 
sic broadenings (which are characterized by the Doppler 
b— parameter in the RATRAN code) take the same value 
everywhere in the numerical model fitting calculations. 

All these averaged spectral profiles present evident red 
profile characteristics, which are indicative of cloud core 
expansions. Our fitting results reveal that L1517B is 
most likely undergoing a core collapse and expansion pro- 
cess as described by the polytropic EECC shock model. 

3.5.3. Several Model Predictions for Cloud Core L1517B 

As predictions for future observations of L1517B, we 
compute spatially resolved relevant molecular line pro- 
files shown in Fig. [5] for molecular line transitions with 
red asymmetries in their central spectra. 

Our model simulation results in Fig. [9] show a gradual 
disappearance of red asymmetries in spectral line pro- 
files and an explicit decrease in overall magnitudes as 
the LOS departs away from the cloud core center. This 
variation trend is attributed to the diminution in both 
cloud density and temperature, as well as optical depth 
(|Gao k. Loull2010T ) . Therefore, it would help to assess the 
influence of optical depths on molecular line emissions by 
making comparisons between model results and observa- 
tions for molecular transitions with spatial resolution. 

4. MILLIMETER AND SUB-MILLIMETER CONTINUA OF 
STAR-FORMING CLOUD CORE L1517B 

While model fittings to observed molecular spectral 
line profiles provide an important means to infer phys- 
ical properties of cloud cores, molecular line emissions 
only present under certain conditions and some molecules 
may be fr ozen onto dust grain surfaces in certain re- 
gions (e.g. iWalmslev et al.M2004l ). Meanwhile, the ra- 
dial profiles of flux intensities at (sub)millimeter wave- 
lengths from dust grain emissions are generally indepen- 
dent of chemical processes, thus can serve as an indepen- 
dent observational diagnostics to examine and constrain 
density and te mperature distributions in star-forming 
clouds (see e.g. iBacmann et al.|[2000l: iShirlev et~al]|2000l: 
lAndre et al.ll2004D . Moreover, as the dust density is al- 
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Fig. 8. — Model fittings of spatially resolved spectral line 
HCO + (l — 0) observations for L1517B. Histograms are Delinglia 
observational data produced by the average of line spectra with 
the same distance from the core center. A telescope efficiency of 
Vmb = 0.62 was used to convert the antenna temperature into 
the main beam brightness temperature T m b. Solid curves present 
LRT fitting simulations averaged over a beam area of 50" X 50" 
based on the same underlying polytropic EECC shock hydrody- 
namic model and molecular abundance pattern. Spectra corre- 
spond to the impact parameter b (distance of LOS from the core 
center) being 0", 20", 50", 71", 100", 112" and 141", respectively. 
The same intrinsic broadening with Doppler b— parameter ~ 0.08 
km s^ 1 and receding velocity of the cloud core « c i ou d ~ 5.73 km 
s _1 as being used previously are adopted in our data fitting pro- 
cedure. 
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Fig. 9. — Red-skewed molecular spectral profiles with spatial res- 
olutions for L1517B. Emissions from HCO+(3 - 2), H 2 CO(2i 2 - 
In), H 2 CO(2n - lio), CS(2 - 1) and CS(3 - 2) are averaged 
over a radius of ~ 20" to simulate the observation of IRAM 30 
m telescope. We present emission spectral profiles with the chosen 
impact parameter b being 0", 40", 60", 72" and 80" respectively to 
reveal the variation of the relative difference in strengths between 
red and blue peaks as well as the overall magnitudes. 
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most independent of depletion holes of molecular abun- 
dances, millimeter and (sub)millimeter continua would 
offer an additional probe to distinguish different models. 
For optically thin dust emissions, the integral form of 
the specific intensity from a spherically symmetric cloud 
globule along a LOS with an impact parameter b is given 

by 

I,{b) = 2^ BAUr)] ^ P ^ 2 dr , (13) 

(e.g. lAdamsl 119911 ) where i? ut is the outer radius of a 
cloud globule, and Td(r), p(r), «v(r) and B„{T) are the 
dust temperature, cloud mass density, specific dust opac- 
ity and the Planck function, respectively. 

4.1. Radio Continuum Emissions at 1.2 mm 

The 1.2 mm wavelength radio continuum mapping of 
L1517B is centrally c oncentrated and ap pears grossly 
spherical (see figure 1 of lTafalla et"aTl20 04). so the spher- 
ical symmetry is a reasonable first-order approximation 
for L1517B. We use the same EECC shock hydrodynamic 
model described in Section 2 to fit the radial profile of 
1.2 mm radio cont in uum emission data ( see figure 2 in 
iTafalla et ail 120041) . iTafalla et all (|2004l) acquired the 
1.2 mm radio continuum mapping from the IRAM 30 
m telescope with a beam size of ~ 11". Therefore, we 
integrate the intensities from model calculations within 
this beam size. Parameters of the EECC shock model 
are shown in Tables 1 and 2, and our theoretical model 
fitting to the observed radial profile of radio continuum 
is displayed in Fig. [10] In the model computations, the 
dust temperature is assumed the same as the gas temper- 
ature adopted for previous molecular spectral line profile 
calculations shown in the Section 4. We utilize the dust 
opacity model proposed by Osscnko pf fc HennineJ (|1994j ) 
with emissions and absorptions included but with scat- 
ters ignored. A typical gas to dust mass ratio of 100 to 1 
is assumed in our cloud model calculations as usual. Nu- 
merical calculations are carried out using the RATRAN 
code with the spherical cloud divided into 256 spherical 
uneven shells. We have also doubled the number of shells 
(i.e. 512) in our calculations (see the dashed curve in the 
top panel of Fig. ITOl) and find that, except for a tiny in- 
crement at small radii, the result generally coincides with 
that produced by 256 shells. Therefore, our adoption of 
256 shells provides sufficient computational accuracy for 
modeling the continuum emissions. Though the density 
and temperature of the cloud increase rapidly near the 
center, they do not lead to a noticeable central peak or a 
prompt increase in the dust continuum since our calcula- 
tion results have been convolved with a telescope beam 
size of ~ 11" which smooths the singularities of density 
and temperature profiles at very small radii. 

4.2. Continuum Emissions for 850 /im and 450 /im 

High-quality sub-millimeter continuum observations 
at wavelengths of 850 /im and 450 /im reveal more 
structural complexities of pre-protostellar cores (e.g . 
Holland et aH[l99l iShirlev et alll2000T : iDve et al J 12003 ). 
As a consistent check, we compute radial profiles for 850 
/im and 450 /im flux intensities using the same EECC 
model and compare them with the observational data of 
SCUBA by Kirk et al. (2005). The same opacity model 



for n v (Oss enkopf fc Henning]ll994l) and the gas-to-dust 
mass ratio of 100 to 1 are adopted in these model calcu- 
lations. 

In order to be compatible with the JCMT beam 
FWHM of ~ 14.8 arcsec at wavelengths 850 /zm and 
^ 8.2 arcsec at 450 /im, our numerical calculations are 
convolved with a telescope beam size of ~ 15" x 15" 
for 850 /im and of ~ 8" x 8" for 450 /im, respectively. 
The fitting result of 850 /im normalized flux density is 
shown in Fig. [10] and the model calculation of the 450 
/im normalized flux density radial profile is shown in Fig. 
[TT1 We find a reasonable agreement of the model results 
with observations of 850 /im continuum inside 10 AU; 
outside this radius, the quality of data becomes poor 
with considerable error bars. The 450 /im flux density 
profile is grossly consistent in intensity with the obser- 
vational data (see figure 3 of Kirk et al. 2005) at cor- 
responding radius, with several substructures being ig- 
nored. As we have not included the influence of the cen- 
tral point source in our calculations, peak flux intensi- 
ties of both 850 /im (~ 155 mJy beam -1 ) and 450 /mi 
(~ 790 mJy beam -1 ) are somewhat lower than those 
observed, namely S^Q k — 170 ± 12 mJy beam -1 and 

Sf 5 ° ak = 920 ± 120 mJy beam -1 (see table 1 of Kirk et al. 
2005). Since the radio observations of 1.2 mm continuum 
and 850 /im data were carried out by two different tele- 
scopes (i.e. the IRAM for 1.2 mm and the SCUBA for 
850 /im) and our model computational result agrees well 
with the 1.2 mm continuum of dust emissions (see Fig. 
[TO]) , the systematic difference of our fitting for the 850 
/zm continuum profile may also be related to differences 
in calibrations between the two telescopes, which could 
be checked by future direct radio observations of 850 /im 
continuum using the IRAM 30m telescope. 

4.3. Column Number Densities and Dust Extinctions 

The near-infrared dust extinction provides an indepen- 
dent diagnostics to trace the column density in molecu- 
lar clouds without involving the dust temperature dis- 
tributions (|Alves. Lada fc Ladal[200l . It serves as an 
important constraint to the density radial distribution 
of our general polytropic EECC shock hydrodynamic 
model, as a complementary yet independent diagnostics 
to various molecular line pro files and dust rad io contin- 
uum emission measurements (ILou fc Gaoll2011lf rB l . More- 
over, column densities in the low-density outer region 
of molecular cloud, where (sub)millimeter dust contin- 
uum emissions may not be easy to detect, can still be 
traced by near-infra red dust extinction measurements 
(|Kandori et al.l 120051 ). As an essential and necessary 
check, we thus compute the column density radial pro- 
file of H2 molecules by simply integrating volume densi- 
ties from the same EECC shock hydrodynamic model. 
Since the column density radial profile is sensitive to 
the outer radius i? ou t of a cloud globule, our model 
fitting to the observed radial profile would identify a 
proper outer radius of the cloud core at i? ut = 2.5 x 
10 4 AU. We then invoke the empirical con version relation 
N(R 2 )/A V = 9.4 x 10 20 cm -2 mag -1 (e.g. iKandori et al.l 

16 Because the temperature profile is not isothermal (as often 
assumed in the literature) in our EECC model, the column den- 
sity radial distribution cannot be directly mapped out by dust 
(sub) millimeter continuum emissions. 
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Fig. 10. — Radial profiles of the 1.2 mm (top) and 850 fim (bot- 
tom) radio continuum emissions from L1517B. For the 1.2 mm 
continuum, the intensity scale is in mjy per 11" beam size and 
the abscissa is the angular distance from the cloud core center i n 
arcsec. Solid squares are data from fig ure 2 of Tafalla etall [I2004T ). 
The solid and dashed curves present the simulation fitting results 
based on the density and temperature profiles of our EECC shock 
hydrodynamic model by dividing the cloud into 256 and 512 shells, 
respectively. For 850 fim, we show the normalized radial profile for 
the flux intensity (per 15" beam size) with the peak flux density 
Sg^g k = 155 mjy beam -1 . Data points from SCUBA observa- 
tions (Kirk et al. 2005) are shown here in a logarithmic scale at 
half beam spacings with lex error bars (~ 17 mjy beam -1 ). The 
abscissa stands for radius from the center in unit of AU on a loga- 
rithmic scale. The scatter of data points, in both the 1.2 mm and 
850 fim continua, at the most outer part of the core (R > 80", 
where the inconsistency between our simulations and observations 
eme rge) is highly likely due to the asymmetry of the core (see figure 
1 of Tafalla et al. 2004). No conspicuous shock discontinuities are 
seen for either 1.2 mm or 850/xm continua, since this shock locates 
at the outer envelope (~ 140 AU) where both the temperature and 
density change slightly due to the shock wave. 

l200l IKirk. Ward-Thompson fc Andrei 12001 . in which 
-/V(H2) is the H2 molecule column number density and 
Ay is the visual extinction, to predict the column den- 
sity and visual extinction for cloud core L1517B using 
our EECC shock model which can reasonably fit other 
currently available observational data. 

In the column density calculation, we convolve numer- 
ical integration results with a typical beam resolution 
of ~ 15". As a comparison and prediction, results of 
both the H2 molecule column density and the visual ex- 
tinction Ay derived from the EECC shock model are 
shown in the lower panel of Fig. [TTJ For example, the 
value of central column density predicted by our EECC 
shock model (~ 2.2 x 10 22 cm~ 2 ) is close to yet a bit 
lower than the model result of Kirk et al. (2005) [i.e. 
A^(H2) C = 4 x 10 22 cm' 2 with a typical error range of 
±20 ~ 30 percent (see table 4 of Kirk et al. 2005). Their 
central "column density" , however, instead of being di- 
rectly obtained from infrared dust extinction measure- 
ments, is actually derived from the 850 fxia radio con- 
tinuum profile modelled by a static isothermal (T = 10 
K ) Bonnor-Ebert sphere and a constant dust opacity 
Kgso = 0.01 cm 2 g _1 . However, with those estimated 
parameters, such an isothermal Bonnor-Ebert sphere for 
L1517B turns out to be dynamically unstable. Besides, 



Fig. 11. — Radial profile for the 450 fim radio continuum emis- 
sions (top); and column densities for hydrogen molecules A^H^) 
and corresponding magnitudes of visual extinction Ay (bottom) 
predicted for observations of L1517B using the EECC shock model. 
For the 450 fim continuum, here we show the normalized radial 
profile for the flux density (per 8" beam size) with the peak flux 

intensity S^g k = 790 mjy beam -1 . Results are shown in a log- 
log scale with the abscissa being the radial distance in AU. For the 
column density and visual extinction profiles, the ordinate on the 
left axis marks the column number density distribution JV(H2) in 
unit of cm - 2 and the ordinate on the right axis indicates the mag- 
nitudes of visual extinction Ay . The abscissa is the radial angular 
distance from the cloud core center in arcsec. These results are 
convolved with a typical observational resolution of ~ 15". 

their central "column density" is model dependent with- 
out being constrained by other available observations of 
L1517B, so we propose here to test and examine their 
prediction as well as ours for the column density radial 
profile of L1517B by dust extinction observations. 

5. CONCLUSIONS AND DISCUSSION 

In this paper, we have systematically investigated the 
physical and chemical (abundance) properties of the star- 
less cloud core L1517B from three complementary ob- 
servational aspects, viz., molecular spectral line profiles 
(some of them are spatially resolved), (sub)millimeter 
continuum emissions and the column density radial pro- 
file through near-infrared measurements of dust extinc- 
tion. We invoke a self-similar polytropic EECC shock 
model with six independent model parameters involved, 
viz. the upstream sound parameter fci, time scale t, poly- 
tropic index 7, coefficients of asymptotic solution condi- 
tions A and B, and the shock location x±, as summarized 
in Table 1 to present reasonable data fittings to various 
asymmetric 'red profiles' together with optically thin sin- 
gle peak spectral profiles, as well as to 1.2 mm and 850 
jiva conti nuum radial profiles from observations (Kirk et 
al. 2005; lTafalla et alj|200i [2001) . In addition, we have 
recently completed a 5 x 5 (with a 50" spacing) spatially 
resolved spectral profile observation of the HCO + (l — 0) 
line transition using the PMO 13.7 m millimeter-wave 
radio telescope at Delingha and made a comprehensive 
comparison with our model results. Furthermore, other 
predicted molecular spectral line profiles with higher spa- 
tial resolutions are expected to be tested by future high 
spatial resolution spectral observations of L1517B. 
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The underlying general p olytropic self-simil ar EECC 
shock hydrodynamic model (Wa ng fc Loul[2008l ) provides 
radial velocity, mass density and temperature structures 
of the cloud globule consistently. Based on the inferred 
EECC model with a set of fitting parameters, we pro- 
vide physical parameters for L1517B. (1) The L1517B 
molecular cloud core is likely evolving in an early phase 
of pre-protostar formation with a small central point 
mass of Mo ~ 0.034 M Q (incapable of thermal nuclear 
burnings) and a current central mass accretion rate of 
M ~ 4.54 x 10" 8 M Q yr" 1 . (2) The total mass of 
the cloud core is estimated as M to t ~ 3.89 M© reveal- 
ing that L1517B b e longs to a low-mass clou d core (e.g. 
ITafalla et all [200l iMcKee fc Ostriken l200l . (3) The 
cloud core L1517B has an expanding envelope with a 
typical outgoing speed of it oxp ~ 0.1 km s _1 and a core 
collapsing at u- ln f ~ 0.2 km s _1 03 (4) There is also an 
expanding shock at radius i? s h = 1-98 x 10 4 AU with an 
outgoing radial speed of ~ 0.25 km s -1 . This outgoing 
shock is initiated when an expanding outflow runs into 
a slowly infalling gas in our dynamic EECC model sce- 
nario. More importantly, we demonstrate that the poly- 
tropic EECC shock dynamic phase of L1517B with a set 
of sensible parameters can give rise to the manifestation 
of various observed molecular spectral line profile s. 

In contrast to the results of ITafalla etaD (|2004L 120061 ) 
and Kirk et al. (2005) (see section 2.1 for details), 
our general polytropic EECC shock hydrodynamic model 
simply depends on physical properties of L1517B in- 
side the observed outer boundary of ~ 150" and self- 
consistently produces thermal as well as dynamic prop- 
erties of the cloud core in the framework of self-similar 
hydrodynamics to fit the molecular spectral line profiles, 
(sub)millimeter continuum radial profiles and to predict 
the dust extinction property and the column density ra- 
dial profile in the plane of sky. 

The uniqueness of our model fitting is ge nerally sup- 
ported by the theoretical work done by IGao fc Loul 
(2010) and Lou & Gao (2011), which discussed the rela- 
tionship between different sets of model parameters and 
the corresponding molecular line spectral profiles. In- 
stead of other possibilities that might cause red-skewed 
molecul ar spectral l i ne pro files as noted in the introduc- 
tion of IGao fc Lov] (l2010h . we conclude that the 'red' 
asymmetries present in several molecular central spec- 
tra of L1517B are most plausibly produced by the co- 
existence of core collapse and envelope expansion mo- 
tions in the cloud. Other scenarios, such as rotation 
an d pure contraction w ithout expansion were suggested 
by ITafalla et al.l ([2004). However, as discussed in Sec- 
tion 13. 5.11 spatially-resolved molecular spectra around 
the core center reveal that this asymmetry is most likely 
an intrinsic property of the cloud dynamics, rather than 
being caused by a rotation. Besides, as the temperature 



profile in our model increases towards the center (a typi- 
cal characteristics for most molecular clouds), it is impos- 
sible to produce molecular line profile s of red asymmetr y 
from a pure contr action (see eq. B4 of|Gao fc Lovfe OlO'). 

As indicated by IBroderick et al.l (|2002t ). there exist ra- 
dial breathing modes generated from nonlinear evolution 
of acoustic radial pulsations on large spatial and tempo- 
ral scales, with long crossing times in protostar-forming 
molecular cloud cores. We advanced the following sce- 
nario for the possible emergence of the EECC phase of 
L1517B. The molecular cloud is undergoing an envelope 
expanding phase of the breathing mode when the cen- 
tral contracting region collapses due to nonlinear insta- 
bilities. This central infalling dynamics is promoted by 
self-gravity and the collapsed regio n expands t owards the 
envelope, similar to the EWCS in iShul (|1977| ). Eventu- 
ally, all gas materials would continue to fall towards the 
center to form a proto-stellar core. This scenario is con- 
sistent with both theo retical (e.g. | Stahler fc Yenll2bl0f) 
and observational (e.g. iLee fc Mversll2011|) pictures con- 
cerning the evolution of internal motions in starless cores. 

We find in our model analysis that a comparison be- 
tween optically thin and thick molecular spectral line 
profiles (or between different energy level transitions of 
the same molecule) serves as an effective and sensible 
method to study the influence of optical depth varia- 
tions. In addition, analyses of the same molecular line 
transition along LOS with distinct impact parameter 
b from the center may further provide an effective ap- 
proach to examine effects of both optical depths and ra- 
dial variations of cloud physical properties (see Section 
3.5). We have prescribed the abundance patterns of dif- 
ferent molecules with central depletion holes, which has 
been widely invoked for molecular line profi les from star- 
forming clouds as radiative diagno sis (e.g. ITafalla et al.l 
l200l [200l iWalmslev etafl l200l . Intensities of 'red' 
peaks appear to be much more sensitive to the size of a 
depletion hole, as compared to 'blue' peaks. This may 
allow us to assess molecular depletions by comparing the 
relative differences between intensities of 'red' and 'blue' 
peaks in their molecular spectral line profiles. 
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17 We show the velocity profile from our model in Fig. [T] these 
typical values only imply the characteristic magnitudes of expand- 



ing and collapsing flow velocities, rather than implying the outgo- 
ing/infalling region has an uniform velocity of this typical value. 
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APPENDIX 

GENERAL POLYTROPIC SELF-SIMILAR HYDRODYNAMIC MODEL 

In spherical polar coordinates (r, 9, </>), general polytropic nonlinear hydrodynamic partial differential equations 
(PDEs) for a spherically symmetric molecular cloud under the self-gravity and gas pressure force are given by 



d P , 1 9 < 2 \ n 



dt ' r 2 dr 
DM ^dM _ Q 
dt dr 

( du du 



dM 2 
— — = Airr p 
or 

dp GMp 
dr r 2 



tt + "7T In — = 
dt dr J\ p-y 



(Al) 
(A2) 
(A3) 
(A4) 



where the mass density p, the radial bulk flow velocity u, the thermal gas pressure p and the enclosed mass M within 
radius r at time t depend on r and t in general; G = 6.67 x 10 -8 dyne cm 2 g -2 is the gravitational constant and 7 is 



the polytropic index. PDEs (|A1[) and (| A2|) are the two complementary forms of the mass conservation. PDE (IA3|) is 
the radial momentum conservation of a molecular cloud unde r the ga s pr essure force and the self-gravity but in the 
absence of random mag netic fields (|Wang fc Loul [20071 l200l . PDE (TMJ requires the specific entropy conservation 
along streamlines corresponding to a general polytropic EoS p — K(r, t)p J with K(r, t) formally related to the specific 
entropy that varies in both time t and radius r in general. 

To derive an important subset of nonlinear self-similar solutions, these hydrodynamic PDEs can be cast into a set 
of coupled nonlinear ordinary differential equations (ODEs) with the following transformation (Wang & Lou 2008), 



= k x ' 2 e 



(the independent self —similar variable x defined) 



a(x) 
AirGt 2 



kt 2 



l p{x) 



4:TTG 



M = 



k 3 / 2 t 3n - 2 m(x) 
(3n - 2)G 



(A5) 
(A6) 



where x is the independent self-similar variable combining r and t in a proper manner, and v{x), a(x), f3(x) and 
m(x) are the dimensionless reduced radial flow speed, mass density, thermal pressure and enclosed mass, respectively. 
Two constants k and n are the sound parameter and scaling index which consistently make x, v(x), a(x), /3(x) and 
m{x) dimensionless and control the spatial and temporal scalings of physical variables in a hydrodynamic cloud. The 
thermal gas temperature is given by the ideal gas law 



T = 



pyrmP 
k B p 



limn 



kt^-^aixy^mixY 



(A7) 



where ks, toh and p are Boltzmann's constant, the hydrogen atomi c mass and the mean molecular weight, respectively. 
For a finite dm(x)/dx as x — > + in the central free-fall solution ('e.g. lWang fc Loull2008f) . the last expression in equation 
16P gives the central mass accretion rate as 



M = fc 3/2 t 3( "~ 1) TO /G 



(A8) 



with 77i being the reduced central enclosed point mass. For n smaller or larger than 1, this central mass accretion 
rate decreases or increases with increasing t, respectively. For n = 1, the central mass accretion rate is constant. 
Substituting self-similar transformation (|A5|) — (|A6I) into PDEs (IA1|) — (|A4|1 . we obtain two coupled nonlinear ODEs 

(A9) 

(A10) 



{nx — v)a — av' — — 2(x — v)a/x , 
f3' /a — (nx — v)v' = (1 — n)v — (nx — v)a/(3n - 



2), 



with /3 = o?m q for 7 7^ 4/3, where parameter q = 2(n + 7 — 2)/(3n — 2). These ODEs are related to the mass 
conservation, the radial momentum conservation and the specific entropy conservation along streamlines, respectively. 
In the limit of x — > +00, we have the asymptotic self-similar solution to the leading orders as 



= Ax~ 2/n + --- , 



v = Bx 1 - 1 ^- 



1 

1-- 

n 



1 A 



(in - 2) 



+ (2n-4)n q - 1 A 1 - n+3nq / 2 



,1-2/r, 



(All) 



(see Lou & Shi 2011 in preparation for more specific comments on the B 2 term) where A and B are two integration 
constants, referred to as the mass and velocity parameters, respectively. 
In the other limit of x — > + and to the leading order, we have the asymptotic central free-fall solution 



a(x) 



(3n — 2)mo 
2^ 



1/2 



v(x) 



2mo 
(3n- 2)x_ 



1/2 



(A12) 
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where the constant mo is the reduced enclosed point mass, representing the dimensionless proto-stellar mass. 

By solving coupled nonlinear ODEs (IA9|) and (IA10[) with analytic asymptotic solutions (jAllI) and (|A12[) as "boundary 
conditions" , and by taking proper care of the sonic critical curve, we can derive the ra dial profiles of velocity, density 
and thermal temperature simultaneously and self-consistently [see IWang fc Loul pOOfih for more details! . 

HYDRODYNAMIC SHOCK CONDITIONS IN THE SELF-SIMILAR FORM 

Self-similar shocks may appear in dynamic molecular clouds and can be constructed within the framework of a 
general polytropic hydrodynamic model. For a hydrodynamic shock, we apply the three conservation laws in the shock 
comoving framework, namely conservations of mass, radial momentum and energy, across the shock front 



01? = o 



[p + p{u s - uf]\ = 



p{u s - u) 3 ~/p{u s - u ) 
2 (7-1) 



-1 2 



= 



(Bl) 



where u, u s , p and p represent the flow velocity, shock front speed, gas mass density and thermal gas pressure, 
respectively. We use a pair of square brackets outside each expression embraced to denote the difference between 
the upstream (marked by subscr ipt '1') and downstream (marked by subscript '2') quantities, as has been done 
conventionally for shock analyses (jLandau fc Lifshitzlll959l ). 

As the parameter k in the self-similarity transformation equation (|A5I) is related to the sound speed which are gener- 
ally different in the upstream and downstream sides of a shock, there are two parallel sets of self-similar transformation 
in the two flow regions across a shock. We set hi — \ 2 k\ with the dimensionless ratio A representing this difference. 
The relatio n x\ = Xx 2 is then required for consistency. With this similar ity scaling relation, hydrodynamic shock jump 
conditions (|B1I) can be readily cast into the following self-similar form ([Wang fc Loull20CM) . namely 



2-n+3nq/2 2q 



a\(nxi —vx) = \a 2 (nx 2 — v 2 ) , 

\g I \2 \2r 2 — n+3nq/2 2q 

vi) q + a.\(nxx — Vi) — A [a 2 x 2 



(nx 2 - v 2 ) q + a 2 (nx 2 - V2) 2 ] 



(nxi — vi) 2 + 



2<l 



(7-1)' 



1 — n+3ng/2 



nxi — v\ 



(nx 2 - v 2 ) 2 + 



27^ 
(7-1)' 



1 — n J r 3nq/2 



(nx 2 —V2) q 



(B2) 
(B3) 

(B4) 



These three self-s imilar shock conditi ons (|B2[) — (|B4[) can be solved explicitly and more relevant details can be found 
in Appendix D of I Wang fc Loul (|2008t ) in the absence of random magnetic fields. 
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